Singularity free cosmological solutions of 
Einstein- Maxwell equations 

Stoytcho S. Yazadjiev*, Ventseslav A. Rizov t 

Department of Theoretical Physics, Faculty of Physics, Sofia University, 
5 James Bourchier Boulevard, Sofia 1164, Bulgaria 

o 
o 

< 



in 



Abstract 

We report on a new two-parameter class of cosmological solutions to the 
Einstein-Maxwell equations. The solutions have everywhere regular curvature 
invariants. We prove that the solutions are geodesically complete and globally 
^ . hyperbolic. 

(N 

o 1 Introduction 

O 

The reasons to study inhomogeneous cosmologies have deep observational and theoreti- 
cal bases. As it is well known the present universe is not exactly spacially homogeneous 
! even at large scales. Moreover, there are no reasons to assume that the regular ex- 

pansion is suitable for a description of the early universe. This strongly motivates 
the study of inhomogeneous cosmological models. They allows one to investigate a 
^ number of long-standing questions regarding the structure formation, the occurrence 

of singularities, the behavior of the solutions in the vicinity of a singularity and the 
possibility for our universe to arise from generic initial data. In particular, on the 
base of some exact solutions it was shown that the nonlinear inhomogeneities could 
regularize the initial singularity giving rise to completely regular cosmologies both in 



general relativity [H]-|15l and in alternative gravitational theories []16|]-|]21|. From a 
purely theoretical point of view, the investigation of nonsingular cosmological models 
gives invaluable insight into the spacetime structure, the inherent nonlinear character 
of gravity and its interaction with matter fields. As a byproduct it also deepens our 
understanding of the singularity theorems, in particular the assumptions lying in their 



base UToH . 



In the present work we present a new two-parameter class of exact singularity free 
solutions of the Einstein-Maxwell equations. 
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2 The exact solutions 



We consider Einstein-Maxwell gravity described by the equations: 

V M F^ = 0, (1) 

We assume that the spacetime admits two spacelike commuting Killing vectors 
<9/<90 and d/dz, which are mutually- and hypersurface-orthogonal. The metric then 
can be written in the well known form 112 211: 



ds 2 = e 2 ^ r \-dt 2 + dr 2 ) + e 2 ^p 2 (t, rW + e~ 2 ^dz 2 . (2) 
We have found the following two-parameter class of exact solutions: 



e 



V(t,r) _ A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at) 



2/ \ .a./r.j.N ' 



cosh (ar) cosh(2at) 



p(t, r) = — cosh(ar) sinh(ar) cosh(2ai), 
a 



= cosh 2 ( 



ar 



A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at) , (3) 



F 03 = 4aAv / T 3 AV 2,/ ' (t ' r) tanh(2at), 



F13 = 4aAv / T 3 A2e- 2V ' ( *' r) tanh(ar), 

where a > and < A < 1 are free parameters. The limiting cases A = and A = 1 
correspond to solutions of the vacuum Einstein equations. The limit a — > gives the 
Minkowski spacetime. The range of the coordinates is: 

- oo < t, z < oo , < r < oo, < (f) < (4) 
The spacetimes described by the solutions (||) have well defined axis of symmetry 



22] . Therefore, it can be said that these spacetimes admit cylindrical symmetry. 

It should be noted that new solutions with the same spacetime geometry and dif- 
ferent Maxwell tensor can be obtained by a duality rotation 

-> cos(#) + *F„ V sin(0), (5) 
where 6 is a constant parameter and * is the Hodge dual. 
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3 Curvature invariants 

The components of the Weyl tensor in the standard null tetrad are: 



3a 2 
cosh 2 (ar) 



2a 2 



8a 2 A 2 
S(t,r) 



+ 16a 2 



A 2 
S(t,r) 
A 2 
' S(t,r) 



1 - 



A 2 



S(t,r) 



cosh 2 (2at) cosh 2 (ar) 
tanh 2 (2at) + tanh 2 (ar) 

2 

[tanh(2at) + tanh(ar)] 2 



1 - 



A 2 



S(t,r) 



[tanh(2at) + tanh(ar)] tanh(2at) 



*i = 0, 



3e 27( *' r % 2 



3a 2 



cosh (ar) 



+ 2a 2 



1 - 



A 2 



8a 2 A 2 



S(t,r) 
A 2 



S(t,r) [ S(t,r)\ 



+16a 2 



+4a 2 



1 - 



A 2 



A 2 



S(t,r) 



cosh 2 (ar) cosh 2 (2at) 
tanh 2 (ar) — tanh 2 (2at) 

tanh 2 (2at) — tanh 2 (ar) 



S(t,r) 



3tanh 2 (ar) - 2tanh 2 (2at) 



^3 = 0, 



3a 2 



where S(t, r) 



cosh (ar) 



2a 2 



1 



A 2 



8a 2 A 2 



S(t,r) 
+16a 2 



1 



S(t,r) 
A 2 
' S(t,r) 

A 2 



+ 



cosh 2 (2at) cosh 2 (ar) 
tanh 2 (2at) + tanh 2 (ar) 



1 2 



[tanh(2at) — tanh(ar)]' 



-8a 2 



1 - 



A 2 



S(t,r) 

A 2 + (1 - A 2 ) cosh 4 (ar) sinh 2 (2at). 



S(t,r) 

[tanh(2at) — tanh(ar)] tanh(2at) 



Further, the components of the Maxwell tensor in the standard null tetrad read 



S(t,r 



-2a\\/l — A 2 cosh(2at) [tanh(2at) — tanh(ar)], 



(13) 



$i = 0, 



(14) 



$ 2 = _ -2aAVl - A 2 ) cosh(2at) [tanh(2at) + tanh(ar)]. (15) 

o [t, r) 

All components of the Weyl tensor and the Maxwell tensor are regular everywhere. 
Since the scalar invariants that can be formed with the metric and the Riemann curva- 
ture tensor are polynomials of these components, all curvature invariants are regular 
everywhere. From the explicit form of the Weyl tensor components one can check that 
the spacetimes are of Petrov type I except at the axis where they are of type D. 



4 Geodesic completeness 

In order to demonstrate the geodesic completeness of the above solutions we have 
to show that all causal geodesies can be extended to arbitrary values of the affine 
parameter. Since the metric functions are even in the time variable we shall investigate 
only future-directed geodesies. Below we consider < A < 1. 

The existence of isometries gives rise to two constants of motion along the geodesies: 



p 



P'(t,r)-, 

-2^(t,r)^£ 

dr' 



where we have denoted by r the affine parameter along the geodesies. 
The affinely parameterized geodesies satisfy 



'dT 

dr 



dr" 
dr 



p 2 (t,r) 



-2ip(t,r) _ p2 e 2^(t,r) = £ 



(16) 



(17) 



where e = 1 and e = for timelike and null geodesies, respectively. Writing d(f>/dr 
and dz/dr as functions of L and P, the geodesic equations for t and r can be written 



in the following form ||14j| , [fl5 



d 


( p 2 7 (t,r) dt 


dr 


{ dr 


d I 


e 2 7 (i,r)fM 


dr~ [ 





r 27(*.r)M(t,r)fikM(t,r), 
,-*r&r)M(t,r)d r M(t,r), 



(18) 
(19) 



where the function M(t, r) is defined by 



M(t,r) 



e + 



p 2 (t,r 



_ e -2^(t,r) + p2 e 2^(t,r) 



1/2 



(20) 



First we consider the null geodesies with L = P = 0. For this case we have 
dt/dr= \dr/dr\ and 



d_ 

dr 



3 2 7 (t,r) 



dt S 

dr 



0. 



After integration, we obtain 



dr 



(21) 



(22) 



where C > is a constant. One can easily see from the explicit form of e 7 ^ ,r - 1 that 
e -2-y(t,r) < ^ As a consequence one finds 



dt ,dr 



< C 



(23) 



dr dr 

and therefore the null geodesies with L = P = are complete. 
Let us turn to the general case when at least one of the constants e, L or P is 
different from zero. In this case we parameterize dt/dr and dr/dr by writing flT^ , 



dt 

dr 
dr 

J? 



-^ {t ' r) M(t,r)cosh(v), 
e -Mt,r)M(t,r)fanh.(v). 



(24) 



Substituting these expressions into the equations for dt/dr and dr/dr, we obtain 
the following equation for v: 



dv 
dr 



-27(*,r-) 



[d t M(t, r) smh{v) + d r M(t, r) cosh(t;)] . 



The above equation can be written in a more explicit form 

dv 1 



dr M(t,r) 



{ J+(t, r) cosh(-u) + J-(t } r) sinh(t;)}, 



where 



(25) 



(26) 



J + (t,r) = 2aetanh(ar) 



1 + 



2(1 - A 2 ) cosh 4 (ar) cosh 2 (2at) 
A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at)_ 
a 3 L 2 [3 tanh 2 (ar) - 1} 
+ tanh 3 (ar)[A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at)] 
+8aP 2 (l - A 2 )[A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at)] tanh(ar), 



(27) 
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X 



(1 - A 2 )cosh 4 (ar) cosh 2 (2at) 
A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at) 
2 [A 2 + (1 - A 2 ) cosh 4 (ar) cosh 2 (2at)] 
cosh 4 (ar) cosh 2 (2at) 
3(1 - A 2 ) cosh 4 (ar) cosh 2 (2at) - A 2 1 tanh(2at). 



J-{t,r) = 4ae 



+2aP 



(28) 



In order for the geodesies to be complete, the functions dt/dr and dr/dr have to 
remain finite for finite values of the affine parameter. In fact, it is sufficient to consider 
only dt/dr since dt/dr and dr/dr are related via (O) and dr/dr cannot become singular 
if dt/dr is not singular. The derivatives d^/dr and dz/dr are regular functions of t 
and r and cannot become singular if t(r) and r(r) are not singular. The only problem 
we could have is when r(r) approaches the value r = for L ^ 0. We shall show, 
however, that r(r) cannot become zero for L ^ 0. 

First we consider the geodesies with increasing r (i.e., v > 0). From the explicit 
form of the functions e 27 ^ and M(t, r) it is not difficult to see that 



~ 2 ^M(t,r) < 



a 2 L 2 



e + P 2 + 



sinh (ar) 



1/2 



(29) 



Therefore, dt/dr could become singular only because of v(t). However, for in- 
creasing r, v(r) cannot diverge since for large t (large r) the derivative dt/dr becomes 
negative as can be seen from (|26|) , (|27|) and (f28|). 

Let us now consider the second case when r(r) decreases (v <0). The geodesies 
with L = reach the axis r = smoothly and then reappear with dr/dr > (v > 0). 
A problem may arise from r = for L ^ 0. However, for L ^ 0, r(r) cannot become 
zero and this can be shown as follows. When r(r) approaches zero the dominant term 
is the one associated with L and the other terms can be ignored. So, for very small r 
the geodesies behave as null geodesies with P = 0: 



^ = e - 27(t ' r) M(r)cosh(i;), 
dr 



dr 
dr 



-2 7 (i,r) 



M(r) sinh(i;) 



(30) 



where 



*d = - e - 2 ^)(9 r M(r)cosh(t;) 

GST 



M(r) = \aL\^ iar) 



sinh (ar) 



(31) 



Hence, one finds 



dr M(r) 



dv d r M{r] 



tanh(f). 



(32) 
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After integration we have 



sinh(ar) = _Dcosh 3 (ar) cosh(-u), (33) 

where D > is a constant. From here one can immediately see that r(r) cannot 
become zero. So, we have proven that our solutions are geodesically complete. In the 
same manner it can be shown that the vacuum solutions for A = and A = 1 are 
geodesically complete, too. 

From the above considerations it follows that every maximally extended null 
geodesic intersects once and only once any of the hypersurfaces t = constant. There- 
fore the hypersurfaces t = constant are global Cauchy surfaces (23| and the spacetimes 
described by the solutions are globally hyperbolic. 

Finally, it is interesting to see which of the assumptions of the singularity theorems 
p4| turn out to be violated. Since the energy and the causal conditions are fulfilled 
it remains to conclude that the spacetimes considered do not contain a closed trapped 
surface. In order to prove this we shall follow the considerations of |2|,[^] and p6| . 
Indeed, let us suppose that the there is one such surface. Then, since the surface is 
compact, there must be a point q where r reaches its maximum. Let us denote it by 
r ma x = R on a constant time hypersurface t = T. For the traces of both null second 
fundamental forms at q, it can be shown that 



K + > V2ae^ R >V 1 + tanh ( 2at ) tanh ( 2ar ) > 0> 

tanh(2ar) 

K- < v ^ Qe -^,T)tanh(2at)t a nh(2ar)-l < Q ^ 

tanni Aav ) 

The traces have opposite sings, and, therefore, there are no closed trapped surfaces. 

In conclusion, we have presented a new (to the best of our knowledge) two- 
parameter class of exact solutions of the Einstein-Maxwell equations. The solutions 
have no curvature singularity. Moreover, they are geodesically complete and glob- 
ally hyperbolic. The solutions can be viewed as explicit examples of how the nonlinear 
inhomogeneities can regularize the singularities yielding completely regular spacetimes. 
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